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Abstract 

We employ the density functional Kohn-Sham method in the local spin-density approximation 
to study the electronic structure and magnetism of quasi one-dimensional periodic arrays of few- 
electron quantum dots. At small values of the lattice constant, the single dots overlap, forming a 
non-magnetic quantum wire with nearly homogenous density. As the confinement perpendicular 
to the wire is increased, i.e. as the wire is squeezed to become more one-dimensional, it undergoes 
a spin-Peierls transition. Magnetism sets in as the quantum dots are placed further apart. It is 
determined by the electronic shell filling of the individual quantum dots. At larger values of the 
lattice constant, the band structure for odd numbers of electrons per dot indicates that the array 
could support spin-polarized transport and therefore act as a spin filter. 

PACS numbers: 73.21.-b, 75.75.+a, 85.35.Be, 85.75.-d 
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I. INTRODUCTION 



Quantum dots or "artificial atoms" , as they are frequently called, confine a few electrons 
on a small conduction electron island, built in (or from) a semiconductor heterostructure. Be- 
ing finite-sized fermion systems, quantum dots can show strong shell effects which determine 
their physical properties. Just like for atoms, quantum dots with closed shells are particularly 
stable, implying "noble" structures for certain numbers of electrons in the dot. Following 
Hund's rules, at half-filling of a shell, orbital degeneracy can lead to spin alignment. This 
was discovered first for small vertical quantum dot samples with circular-parabolic shape by 
Tarucha et al. 1 . The experimental findings were later theoretically confirmed by electronic 
structure calculations using mean field methods as well as quantum Monte Carlo techniques 
or even a numerical diagonalization of the full many-body Hamiltonian (see Reimann and 
Manninen 2 - for a review) . 

Experimentally, few-electron quantum dot structures where the shell effects on magnetism 
could be observed, are challenging to fabricate. One example where spontaneous mag- 
netism has been found, are one-dimensional quantum point contact constrictions formed in 
a gate-patterned heterostructure^^. The intrinsic magnetic properties of these nanostruc- 
tures have drawn much attention recently due to their potential applicability in spintronics 
devices 6 . Quantum point contacts^ and single quantum dots^ were found to have spin filter- 
ing capabilities, with a possibility to serve for either generating or detecting spin-polarized 
currents. 

Arranging many quantum dots in a lattice, one can build artificial crystals with designed 
band structur e) 10 ' 11 which can be manipulated for example by tuning the inter-dot coupling 
and the number of confined electrons in the single quantum dots. The dot lattice does not 
suffer from structural deformations, which has the advantage that it can be designed freely 
without having to consider lattice instabilities 12 . 

Fabrication of a quasi one- dimensional artifical crystal consisting of a sequence of a few 
quantum dots was suggested by Kouwenhoven et a/.— already in 1990. They observed 
oscillations in the conductance as a function of gate voltage, arising from the mini-band 
structure in the periodic crystal. Small dots in well-ordered lattices could be synthesized by 
self-organized growth^. A particularly interesting artificial lattice structure is the Kagome 
lattice, due to the possibility of flat-band ferromagnetis m 12 i 14 i 15 i 16 i 17 . Shiraishi et a/.— have 
pointed at the importance of these structures for fast processing and high-density storage 
of information. 

For square lattices, Koskinen et al. 18 showed within the density-functional scheme that 
few-electron quantum dot lattices have a rich magnetic phase diagram, depending on the 
lattice constant and electron number. Related observations have been made also within the 
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Hubbard mode l 12 ' 19 ' 20 . 

In this paper, we investigate the electronic and magnetic properties of quasi one- 
dimensional quantum dot arrays. We suggest that such linear quantum dot chains could, in 
fact, lead to single-spin conductivity. 

In our model the single quantum dot confinement is provided by a rigid Gaussian-shaped 
background charge distribution. At the single dot centers, this potential is approximately 
parabolic. The band structure and the magnetic properties depend on the lattice constant, 
a, and the number of electrons per dot, N. Here, conductivity of the dot chain is only 
considered by observing whether there is a band gap at the Fermi- level or not, which allows 
a qualitative understanding. 

At small values of the lattice constant, the single dots overlap, forming a non-magnetic 
quantum wire with nearly homogenous density. As the confinement perpendicular to the 
wire is increased, i.e. the wire is squeezed to become more one- dimensional, the ground 
state is a spin density wave caused by a spin-Peierls transitio n 21 ' 22 ' 23 . Magnetism sets on as 
the lattice constant is increased. It is determined by the shell structure of the individual 
dots: the arrays are non-magnetic insulators for closed single-dot shells at N = 2 and 6. At 
half-filled shell (N = 4) the spin of the dot is determined by Hund's rule and the array is an 
antiferromagnetic insulator. Ferromagnetism is observed both at the beginning and the end 
of a shell (here N = 3, 5). The spin-up and spin-down bands are separated by the exchange 
splitting. At sufficiently large lattice constant a one observes a gap between these bands. In 
this case the current would be carried by a single spin only, acting as a spin filter. 



II. THE COMPUTATIONAL METHOD 



In order to model the one-dimensional quantum dot array, we consider interacting elec- 
trons moving in two dimensions in a rigid periodic background charge distribution eps- 
The background charge number per unit cell is chosen to match the electronic charge of 
the unit cell in order to ensure overall charge neutrality. We employ the Kohn-Sham 
method with periodic boundary conditions. The Kohn-Sham orbitals are of Bloch form 
VWr(r) = exp(ik ■ r)w nk(J (r), where n labels the band, a = (|, f ) is the spin index and the 
wave vector k is confined into the first Brillouin zone. The periodic functions u n ko-(r) satisfy 
the Bio ch-Kohn- Sham equations 

h 2 

where the periodic effective potential is 

*/w=/ £ ^^r*' +,4 - WrU(r)1, (2) 
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p is the electron density and £ = (p-f — is the polarization. In the local spin- 

density approximation we use the generalized^ Tanatar-Ceperley2£ parameterization for 
the polarization-dependent exchange-correlation potential v° c [p(r), £(r)]. In the band struc- 
ture calculation, the functions u n y_ a {v) are expanded in a basis with 11 x 11 plane waves. 
For one-dimensional systems, the wave vector reduces to wave number for which we chose 
an equidistant 19-point mesh in the first Brillouin zone. The self-consistent iterations were 
started with anti-ferromagnetic and ferromagnetic initial potentials. Small random pertur- 
bations were added to the initial guesses in order to avoid convergence into saddle points 
of the potential surface. In addition, we use an artificial temperature to allow fractional 
occupation numbers for nearly degenerate states at the Fermi level. We noted that by de- 
creasing the temperature the amplitudes of the spin-density and the average spin per dot 
become somewhat higher for small lattice constants. Nevertheless, we must emphasize that 
the temperature is low enough not to affect the ground-state. The statistical occupations 
merely help occupying degenerate levels to ensure convergence. We use effective atomic units 
with Hartree Ha = m*e 4 /ft 3 (47reoe) 2 for energy and the Bohr radius a* B = h 2 47te e/m*e 2 for 
length, where m* is the effective mass and e the dielectric constant of the semiconductor 
material in question. 

III. MAGNETISM IN A ID QUANTUM DOT ARRAY 

Studying magnetism in a one- dimensional array, the simplest geometry to choose for 
the unit cell is a rectangle with two quantum dots per cell. These dots lie in a row along 
the x axis of the cell, one in the center and one crossing periodically the edge of the cell. 
The confining potential is modeled by a periodic positive background charge distribution 
described by a sum of Gaussians centered at lattice sites R = a(n x , 0), n x = 0, 1, 2, . . ., 



where r = (x, y) is a two-dimensional position vector. A single Gaussian carries positive 
charge Ne with density 1/irr 2 at the center. The parameter r s determines the average 
electron density at the center of the dot. Throughout this paper we use the value r s = 2 a* B 
which is close to the equilibrium density of the two-dimensional electron gas. The bottom 
of the confining potential provided by the background charge distribution is harmonic to a 
good approximation. Since there are two quantum dots in the unit cell, the electronic levels 
are split into bonding and anti-bonding bands. As a consequence, for both spins there are 
two ls-bands, four p-bands, six 2sld-bands and so on. In a one-dimensional quantum dot 
array one can have a smooth transition from the tight-binding description to the nearly-free 
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a=3.14 a=5.0625 a=6.97 a=8.55 a=10.125 



FIG. 1: Lowest bands at selected values of the lattice constant a for a quantum dot array with 
three electrons per quantum dot (in atomic units, see text). The spin-down bands are plotted in 
blue color, and spin-up bands are plotted in red. The dashed green line indicates the Fermi-level 
fixed at zero energy. 

electron picture simply by varying the lattice constant a. 

Figure Q shows the bands for N = 3 with different inter-dot separations. Spin-up and 
spin-down bands are plotted in red and blue color, respectively, and the Fermi-level is fixed 
at zero energy. The spin degeneracy is lifted by the exchange splitting. For very large values 
of the lattice constant a, the electron densities of the single dots hardly overlap, and the 
dots are isolated. 

The energies of ferromagnetic and anti-ferromagnetic solutions are nearly degenerate as 
the local approximation is unable to distinguish between them. Furthermore, the bands 
are flat with band gap energies approximately equal to the single dot level spacings. Even 
though the Fermi energy stays inside a band, the dot array becomes an insulator due to a 
diminished hopping probability between the single dots in this limit. 

By decreasing a, i.e. by bringing the quantum dots closer to one another, the band dis- 
persion increases. The bands corresponding to some specific quantum dot level are bunched 
and the bunches are separated by energy gaps, which is demonstrated in figure [T] for lattice 
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FIG. 2: Magnetism in a linear chain of quantum dots as a function of the number of electrons per 
dot, and the lattice constant. Blue color corresponds to the conducting, yellow to the insulating 
phase. Green indicates the phase where only one spin is conductive. 

constant a = 10.125 a* B . By decreasing a further, the band gaps close. For sufficiently small 
a the quantum dots overlap strongly, which leads to an essentially homogenous quantum 
wire with a Gaussian cross-section. In this nearly-free limit the transverse motion sepa- 
rates from the longitudinal one. Consequently, the transverse states are quantized by the 
Gaussian-shaped well, while the longitudinal states remain "free" with parabolic dispersion. 
This is reflected in the band structure, showing nearly equidistant sub-band parabolae where 
the nth sub-band for a given k corresponds to a Kohn- Sham- Bio ch orbital with n — 1 nodes 
in the transverse direction. We note also that the higher bands have parabolic dispersion 
at longer inter-dot separations than the lower ones due to the longer spatial extent of high- 
energy orbitals. ^From figure Q we note that the second transverse sub-band is occupied 
at a = 5.0625 a* B while at a = 3.14 a* B the Fermi-level reaches the third sub-band. Hav- 
ing now also the higher transverse modes occupied, the quantum dot chain becomes quasi 
one-dimensional . 

Figure 121 shows the magnetism of the quasi-lD quantum dot array as a function of electron 
number per quantum dot and lattice constant a. The colors indicate the regions where the 
array is conducting (blue) or insulating (yellow). The green bars indicate regions where the 
Fermi-level resides solely on a single spin band. The arrows indicate the spin arrangement 
in the array. 

For a single electron per quantum dot, N = 1 only the bonding s-band is filled. Due 
to the exchange splitting of the single dot levels, the bonding and anti-bonding bands are 
separated by an energy gap and the array shows anti-ferromagnetic order. Figure 01 shows 
that the average spin per dot, calculated by integrating the spin density over a single dot, 
drops gradually from 1/2 to as the lattice constant is decreased. The band gap and thus the 
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Lattice constant a 

FIG. 3: Spin per dot for iV=l,3,4 and 5 as a function of lattice constant. 

antiferromagnetism persists down to very small values of the lattice constant. At the closed 
shell N = 2 the bonding and anti-bonding ls-bands are filled leading to a non-magnetic 
insulator. The transition from a tightly bound insulator to a nearly free metal occurs at the 
lattice constant a ~ 6 a* B , when the gap between the Is- and p-bands closes up. 

Next, the p-bands are occupied. There are two degenerate p-orbitals for a single dot 
giving rise to two bonding bands and two anti-bonding bands for both the spins as shown in 
figure |3J The orbitals with density lobes oriented along the wire yield higher dispersion than 
the ones perpendicular to it. For N = 3 there is one p-electron per quantum dot, which 
triggers ferromagnetism. An example of the total electron and spin densities is shown in 
Fig. The levels with majority spin are lower than the ones with minority spin as a result 
of the exchange splitting of the energy bands. The density in the array increases as the dots 
are brought closer. Consequently, the kinetic energy becomes the dominant contribution to 
the total energy. ^From figure Q we note that at a = 5.0625 a* B the dispersion is parabolic 
and the spin degeneracy is restored. However, there is a competition between the kinetic 
and exchange energies. It turned out that at a = 3.14 a* B a small spin-splitting is re-gained. 

Since the Fermi-level is bound to the p-band region, the array with three electrons per dot 
remains conductive at all lattice constants. Figure H] shows that the bands of the minority 
spin are pushed up in energy by exchange splitting. As a consequence, just before the 
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(a) N=3, a=1 3.05 (b) N=5, a=1 5.3 




FIG. 4: lp-bands and lowest 2sld-bands for (a) N = 3 at a = 13.05 a* B and (b) N = 5 at 
a = 15.3 a* B . The blue lines are the spin-down bands, the red lines the spin-up bands. The dashed 
green line indicates the Fermi level. The spin-up and spin-down bands are separated by a gap 
which leads to single-spin conductivity. 




FIG. 5: Total density and spin density for N = 3 at lattice constant 13.05 a* E 
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insulating phase, when the bandwidths are relatively narrow, the Fermi-level for the minority 
spin lies in the band-gap but the majority spin remains conductive. A similar behavior is 
found in the case of N = 5 since there are three p-electrons and the shell is almost filled. 
However, this time the minority spin is conductive because the levels with fewer electrons 
are pushed to higher energies. The spin-dependent conductivity might open an intriguing 
opportunity to use the linear quantum dot chain as a spin filter. 

At half-filled p-shell (N = 4) Hund's rule leads to maximized spin in an isolated quantum 
dot. Indeed, the spin per dot for N = 4 is at its maximum (1.0) at lattice constant a « 18 a* B 
and it decreases gradually with a as the spin densities "spill" into the other dots. Due to 
Hund's rule and the exchange-splitting, the anti-ferromagnetic ordering of spins is favored. 
Since a gap is formed at the Fermi-surface, the array is insulating until a transition to 
a nearly homogenous wire occurs. Finally for N = Q, the p-shell is filled and the linear 
quantum dot chain remains non-magnetic at all values of a. 

IV. SPIN-PEIERLS TRANSITION IN A HOMOGENOUS QUANTUM WIRE 

At small values of the lattice constant a the quantum dots overlap significantly, forming 
a homogenous quantum wire with a Gaussian cross-section. Let us look at this limit more 
closely. Consider a quantum wire with a Gaussian cross-section closed in a rectangular unit 
cell. The background charge distribution is chosen to be 



where r\ is the one-dimensional density parameter. The wire lies along the x-axis and its 
width is measured by the full width at half maximum 2^2 In 2a. Since there is no definite 
lattice parameter for the wire, the length L of the unit cell is chosen such that pB integrates 
to the desired charge Ne, thus we have L = 2r]. D N . We have chosen four electrons in the 
unit cell (N = 4) and fixed r\ D = 2 a* B . In addition, we define parameter Cm = 2r\ D /a 
to describe the ratio of the average inter-electron separation and the width of the wire: 
with increasing Cm the wire becomes narrower. Consequently, the energies of the higher 
transverse modes are pushed up in energy. 

Figure |H1 shows band structures of homogenous quantum wires for selected widths. For 
Cm = 2, the dispersion is parabolic and the Fermi-level lies close to the second transverse 
sub-band. In this case the wire shows no magnetism. Antiferromagnetism sets on at Cm = 4, 
as the spin-Peierls transition occurs. The ground state is a spin density wave with wave 
length of L/2 = r] D N = 8 a* B . The spin-Peierls transition opens a gap at the Fermi- level 
and turns the wire into an insulator. The amplitude of the spin density wave increases with 
increasing Cm- 




(4) 
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FIG. 6: Lowest bands at selected values of width parameter C\d for a quantum wire with four 
electrons per unit cell. The dashed green line indicates the Fermi-level fixed at zero energy. 

V. SUMMARY 

We studied the electronic and magnetic properties of one-dimensional arrays of few- 
electron quantum dots. The spin per dot, and thus the magnetism of the array, depends 
on the shell filling of the individual dots and the inter-dot coupling. Furthermore the band 
structure of chains with open-shell dots suggests that conductivity could become spin de- 
pendent at suitable values of the lattice constant. 
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